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We report on new data from additional zero temperature simulations of QCD with 12 flavors. 
This is a continuation of previous studies using the DBW2 gauge action and naive staggered 
fermions. With the use of the force gradient integrator and a multiple-quark-mass preconditioned 
HMC, we have done simulations with input quark masses from m^ ~ 0.003 to niq = 0.008. We 
have observed a metastable, first order, bulk transition that occurs at small input quark masses. 
As the quark mass increases, this first order bulk transition ends at a second order critical point, 
and, for still heavier quark masses, becomes the cross-over we have previously reported. We 
present measurements of hadron masses, decay constants and other low energy observables in 
the small quark mass region on the weak coupling side of the bulk transition. Our results show 
that the behavior of the system is still consistent with spontaneously broken chiral symmetry. 
We also discuss a preliminary investigation into the behavior of the bulk transition itself. We 
have found that, as the system approaches the second order critical end point, the scalar singlet 
meson becomes lighter Thus it appears that the critical endpoint corresponds to a continuum 
limit theory only involving scalars and, following known triviality arguments, this is likely a free 
field theory. The presence of this critical endpoint could influence scaling of lattice observables 
in the conventional continuum limit. 
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(a) 8 flavors, m^ = 0.015 
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(b) 12 flavors, m„ = 0.01 



Figure 1: The evolution of (v^V'^) with 8 quark flavors (left) and 12 quark flavors (right), containing ordered 
(bottom curves) and disordered (top curves) starts, both with a lattice size of 16^ x 32. 

1. Introduction 

The zero temperature phase of QCD with many degenerate massless quark flavors has been 
a heated topic for a number of years. It is important to establish if an infrared fixed point is 
dynamically generated by a non-abelian gauge theory. It enables all the theoretical possibilities 
with a conformal symmetry (if the fixed point exists) or a walking coupling constant (if the number 
of flavors is less than required for a fixed point), which leads to some exciting theories as candidates 
for beyond the standard model physics. 

Lattice methods are indispensable for determining the critical number of quark flavors required 
for the fixed point, as it is a non-perturbative phenomenon. Various groups are working on the sub- 
ject. As the situation is new and unlike QCD in nature, we need to be careful in our interpretations 
of the numerical results from lattice calculations. 

We have been using the DBW2 gauge and the naive staggered fermion action to simulate the 
SU(3) gauge theory with 8 and 12 flavors. In this report, we will discuss our investigation into a first 
order bulk transition separating strong and weak couplings, for small quark masses, and the second 
order critical point, where one expects the bulk transition to end, as the quark mass is increased. 
We argue that such a critical point (a lattice artifact) can alter the scaling behavior of the hadronic 
observables, if simulations are done in regions of the parameter space where the critical point has 
influence, and may have an impact on the running of the coupling constant directly measured on 
the lattice. 



2. First Order Bulk Transition 

A bulk transition has been found and studied in an earlier paper[l], for lattice QCD with 8 
quark flavors, using the Wilson gauge action and naive staggered fermions, with the inexact R 
algorithm. The bulk transition was interpreted as a lattice artifact. It is independent of the lattice 
volume and temperature, and it is not associated with a deconfinement transition or chiral symmetry 
restoration. Recently, using the exact RHMC algorithm, we have rediscovered the bulk transition 
with the same gauge and fermion action as in the earlier study. We found it at j8 = 4.59 and 
ntq = 0.015 with a lattice size of 16^ x 32, shown in figure la. Such a bulk transition, however, 
did not occur within the parameter space we have surveyed in our recent studies [2, 3, 4] with the 
DBW2 gauge and naive staggered fermion action. 

One might expect the lattice artifact bulk transition to become stronger with 12 quark flavors. 
The increase in the number of flavors makes the renormalized coupling run more slowly, and a 
higher lattice cutoff is required to approximate the continuum. In other words, at a comparable 
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Figure 2: (</</) (left) and /;i (right) with 12 flavors. Dashed lines indicates the meta-stable signal from the 
first order bulk transition, seen for both (5 = 0.46 and 0.4626. 
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(b) j8 = 0.4626, nig = 0.008 



Figure 3: The evolution of (<//</) with DBW2 gauge and naive staggered fermion action. Two combina- 
tions of j3 and nig show meta-stable evolutions with separate curves for ordered-starts (bottom curves) and 
disordered-starts (top curves). The metastability is present for two different volumes. 



lattice spacing, the lattice gauge field becomes rougher with increasing number of flavors. The 
lattice consequently becomes a poorer approximation of the continuum. Observation of a bulk 
transition is mentioned in a study of a lattice Schrodinger functional calculation [5]. For 12 flavor 
QCD, with the same Wilson gauge and naive staggered fermion action, we have located the position 
of a first order bulk transition at j3 =4.1 and niq = 0.01 with a lattice size of 16^ x 32 and the exact 
HMC algorithm, as shown in figure lb. 

In our recent study of QCD with 8 and 12 flavors [3] with the DBW2 gauge and the naive 
staggered fermion action, it is intriguing to see that there is no meta-stability signal, but only a 
quark mass dependent, quick cross-over region. Nevertheless, the strong dependence on the quark 
mass suggests that a first order bulk transition may occur with smaller quark masses than we had 
been able to simulate with. Simulations with smaller quark masses became viable with algorithmic 
improvements. Using the force gradient integrator, which is accurate up to the 4* order, combined 
with multiple Hasenbusch mass preconditioning, we were able to use the exact HMC algorithm to 
explore the smaller quark mass region. 

A first order bulk transition was indeed found with 12 quark flavors at smaller quark masses, 
using the DBW2 gauge and naive staggered fermion action. The location of the bulk transition 
is shown in figure 2. It shows the first order, meta-stable points (those have two points at the 
same j3 and iriq), along with other values' from our earlier simulations. In the evolution of (yy)^ 
shown in figure 3, we have found meta-stable signals at two different couplings and quark mass 



'in this paper, f,i is normalized to be about 93MeV for QCD in physical units. Values of /^ in our previous 
reports [3, 4] are off by a factor of 2, because of a misunderstanding of the convention of the CPS. 
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combinations. The meta-stable evolution lasted for a few thousand trajectories. This first order 
bulk transition has no visible volume dependence, as seen from comparing two different volumes, 
16^ X 32 and 24^ x 32. Therefore, it is unlikely to be related to the thermal transition. In addition, 
as discussed previously [3], there is no visible change in the values of Polyakov loops, and there are 
clearly non-zero string tensions found in the inter-quark potentials at both sides of the transition. 
We can conclude that this bulk transition is not a deconfinement transition. On the other hand, 
the chiral symmetry is spontaneously broken on both sides of the bulk transition, as argued in our 
previous report [3]. The signal from all the hadronic observables is clear and well justified at the 
strong coupling side of the transition. On the weak coupling side, the linearity of m| versus mq 
suggests the existence of a dimensionful quantity B, as in leading order chiral perturbation theory 
where m^ = IBniq. Consequently, with a non-zero f-j^ in the chiral limit, the GMOR relation, 
/^ni| = 2m{\j/Y), predicts a small but non-zero chiral condensate in the chiral limit, which is 
consistent with out data. This argues that the weak coupling side of the bulk transition also has 
the chiral symmetry spontaneously broken. More of the spectrum results will be shown in the next 
section. In conclusion, the bulk transition does not restore chiral symmetry from the data we have 
collected. 

Our simulations suggest the following phase structure. There is a first order bulk transition 
line in the /3 versus niq plane, at smaller masses, riiq < m^ At larger masses, niq > niq, the bulk 
transition becomes a mass dependent rapid cross-over. In addition to the usual Gaussian fixed point 
at j3 = oo and niq = 0, where the lattice gauge field approaches the continuum limit asymptotically, 
the results from our study of the system with 12 quark flavors suggest that there is an additional 
critical point at j3 = j8^ and triq = m^. It is the end point of the first order bulk transition, a 
second order phase transition point, where certain correlation lengths diverge. It also represents a 
continuum theory. 

The position of the critical point depends on the lattice action. In our case, the change from 
the unimproved Wilson gauge to the DBW2 gauge action moves the critical point to smaller quark 
masses in lattice units. 

With 8 flavors, it is also possible that the bulk transition does not go away completely, after 
we have switched to DBW2 from unimproved Wilson gauge. The bulk transition and the second 
order critical point may persist at a smaller quark masses away from the quark mass region we 
have analyzed. This possibility, however, does not invalidate our conclusions about the phase of 
the QCD system with 8 flavors. The spontaneous breaking of chiral symmetry and the asymptotic 
freedom holds at both side of the bulk transition. 

Nevertheless, it is plausible that this additional critical point (different from /3 = oo) is capable 
of affecting the scaling of the hadronic variables and also changing the running of the coupling 
constant directly measured on the lattice. More careful and thorough studies are needed to investi- 
gate the effect on the running of the coupling constant. Yet, continuing our recent studies, we are 
able to identify, detailed in the next section, the effect of this bulk transition critical point in the 
meson spectrum. 

3. The Meson Spectrum and the Scalar Singlet Meson 

Results for some meson masses are shown in table 1 for jS = 0.46 and j8 = 0.4626, respectively, 
at mq < tiiq. Wall2Z sources are used for the propagators, except for /;;;:, which we measured from 
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24-^x32 
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0.1191(14) 
0.2080(28) 


0.24368(28) 
0.21068(14) 


0.565(16) 
0.882(97) 


0.756(12) 
1.151[13] 


0.995 [94] 
1.61[51] 



0.4626 



0.008^ 24^x32 0.1442(15) 0.26973(18) 

0.008^ 24^x32 0.1652(16) 0.25926(16) 

Table 1 : Some of the meson masses measured in 
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Figure 4: Masses of the n, p and oi mesons versus input quark mass, nig, at )3 = 0.46 and j3 = 0.49. 



point to point propagators. In the table, the lable "^ or ^ on the quark mass indicates an ordered 
(weak coupling side) or disordered (strong coupling side) start, respectively. Most of the errors are 
statistical, except those enclosed by angular brackets have been increased to include a systematic 
variation caused by varying the fitting range. 

Figure 4a shows that, on both sides of the bulk transition, there is the profound linearity of 
m| = IBiTiq. This suggests that the chiral symmetry is spontaneously broken on both sides of the 
bulk transition. It is confirmed in the mass splitting between p and a^ mesons, as depicted in 
figure 4b. 

There is still one thing left unexplained. It is well known that the correlation length diverges 
as a system approaches a second order critical point, which means some mass in the system with 
12 flavors should vanish at the critical point. We can see from table 1 that, while the system goes 
from strong to weak coupling through the bulk transition, in lattice units, the masses of a, p and ay 
mesons drop, and the pion mass increases. This means that, going from strong to weak coupling, 
the lattice constant or the scale of the system decreases. However, none of the masses here appears 
to be vanishing toward the critical point, from j3 = 0.46, niq = 0.006, to j3 = 0.4626, niq = 0.008. 

With further study, we have found that it is the scalar singlet meson whose mass appears to 
be vanishing near the critical point. Remember that the a meson masses shown in table 1 are only 
from the connected propagators. To measure the disconnected scalar propagators, we employed 
random Z2 volume sources. We have found marginal advantages with the Z2 against the Ul random 
source. We used 40 random sources for each gauge configurations to estimate the quark loops with 
G(t) = {m(p^(p), as proposed by the authors of the paper [6], where M(p = rj with T] the random 
source and M the fermion matrix. The singlet propagators, then, can be assembled as 
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Figure 5: Effective masses of the scalar singlet meson. Gauge configurations with ordered and disordered 
starts respectively correspond to weak and strong coupling sides of the transition. 
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(a) Ordered-start, 440 configurations in total 



(b) Disordered-start, 320 configurations in total 



Figure 6: Error of the effective masses of the scalar singlet meson versus the jackknife block size at j3 



0.4626 and nig = 0.008. 



where O {t) is the connected point to point propagator, (• • • )t is average over time slices, and (• • • ) 
is ensemble average. In the equation, a factor of Nf is removed from the usual definition, and 
the factor of 4 comes from the staggered Dirac operator. We used a single mass state to fit the 
propagator to obtain the singlet mass, as the opposite parity state is numerically invisible to us. 

Figure 5 shows the effective masses of the scalar singlet meson at the bulk transition. The value 
of nieff at T is obtained from C(t) andC(T-|- 1). It shows a bizarre behavior for all four propagators. 
When the temporal separation becomes larger, the effective mass starts increasing. However, it 
is understandable. The singlet propagator becomes noisier and more difficult to measure with 
increasing temporal separations. The expected behavior is demonstrated with the ordered-start 
configuration at j3 = 0.46 and riiq = 0.006, shown in figure 5a. We can see that although the 
expectation value increases with increasing T, the error increases at the same time. 

The errors of ntes are underestimated because of the long autocorrelation time. This is ex- 
pected as the system approaches the second order critical point. As shown in figure 6, the error bar 
grows as one blocks the data in larger bins. It is important to keep in mind that our errors for the 
scalar single meson may be underestimated. 

The mass of the scalar- singlet meson was measured via fitting the propagators. The table 
at the left side of figure 7 shows the fitting range and the results. The statistical eiTor quoted 
here were estimated using the jackknife method with a block size of 8 configurations. It is very 
exciting to see that the scalar singlet meson becomes really light at the point of the bulk transition 
at /3 = 0.4626 and m^ = 0.008, since it is closer to the second order critical point. If we apply 
the scaling of universality class as m^v ~ C(m^ — m^)'/^, shown in figure 7, we can estimate the 
critical quark mass at the end point to be at m^ = 0.00883(45) from ordered-start configurations, 
and m^ = 0.00869(29) from disordered-start configurations. 
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Figure 7: Masses of the scalar singlet meson at the first order bulk transition, and the extrapolations to the 
critical point. 

4. Conclusion 

We have found a first order bulk transition with 12 light quark flavors, using the DBW2 gauge 
and naive staggered action. All the hadronic observables are consistent with the system having 
spontaneously broken chiral symmetry and asymptotic freedom on both sides of the bulk transition. 

The scalar singlet meson becomes lighter as the system approaches the second order critical 
point (the end point of the first order bulk transition), and is expected to have vanishing mass as 
other masses or hadronic scales remain finite at the critical point. This suggests that the continuum 
limit of the system at this critical point is a theory of scalars and hence one expects it to be a free 
theory, from the general arguments of the triviality of scalar field theories. 

The scaling behavior of the hadronic observables is expected to be altered near this critical 
point, as opposed to the continuum limit governed by the Gaussian fixed point at j8 = oo. Care must 
be exercised to do any scaling test. In particular, measuring the sigma mass while approaching the 
conventional continuum limit should indicate how much influence this critical endpoint has at a 
particular point in parameter space. 

The running of the gauge coupling directly measured on the lattice is expected to be altered 
by this critical point, as it potentially deviates from the behavior in the weak coupling limit. More 
study is required to quantify this effect. 

Our calculations were done on the QCDOC and NY Blue at BNL. X.-Y. Jin wants to express 
the sincere gratitude especially toward all the brilliant minds behind the QCDOC project. This 
research utilized resources at the New York Center for Computational Sciences at Stony Brook 
University/Brookhaven National Laboratory which is supported by the U.S. Department of Energy 
under Contract No. DE-FG02-92ER40699 and by the State of New York. 
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